Mesons Above The Deconfining Transition 
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We analyze temporal and spatial meson correlators in 
quenched lattice QCD at T > 0. Above T c we find different 
masses and (spatial) "screening masses", signals of plasma 
formation, and indication of persisting "mesonic" excitations. 

PACS numbers: 12.38. Gc, 12.38.Mh 

With increasing temperature we expect the physical pic- 
ture of QCD to change according to a phase transition 
where chiral symmetry restoration and deconfincmcnt 
may simultaneously occur. For model independent non- 
perturbative results one attempts lattice Monte Carlo 
studies Since in the Euclidean formulation the 0(4) 
symmetry is broken at T > (see, e.g., ||), physics ap- 
pears different, depending on whether we probe the space 
("er": x) or time ("t": i) direction: the string tension, 
e.g., measured from era Wilson loops does not vanish 
above T c , in contrast to the one measured from <jt loops. 
Therefore we need to investigate hadronic correlators with 
full "space-time" structure, in particular the propagation 
in the Euclidean time. The latter, however, poses special 
problems because of the inherently limited physical length 
of the lattice in the time direction l T = 1/T. We shall 
discuss briefly this question and introduce our procedure. 

1) Lattice problems: Large T can be achieved us- 
ing small N T = l T /a (a: lattice spacing), however this 
leads to systematic errors ||. Moreover, having the t- 
propagators at only a few points makes it difficult to 
characterize the unknown structure in the correspond- 
ing channels. To obtain a fine ^-discretization and thus 
detailed ^-correlators, while avoiding prohibitively large 
lattices (we need large spatial size to avoid finite size ef- 
fects, typically l a ~ 3Z T ), we use different lattice spacings 
in space and in time, a a j a T = £ > 1. For this we employ 
anisotropic Yang-Mills and fermionic actions j| : 

Sym = -(/3/3)( 7 g 1 ReTrD CTCT + 7G ReTrD CTT ) (1) 
S F = {2 Ka )- 1 *WV, k- 1 = 2(m + 3 + 7f), (2) 
W = l-K a { ^ T+U^ + 1F T+U A T 4 + "h.c." ) (3) 



(r± = i± 7Ml7 2 



1 , TV, are lattice translation operators 



and mo the bare quark mass). £ is determined from T ~ 
correlators F ( "calibration" ) by requiring isotropy in the 
physical distance: F (7 (z) = F T (t — £z). In a quenched 
simulation at some 7g, £ is fixed by the Yang-Mills cal- 
ibration, then jp is tuned to give the same £ for T ~ 
hadron correlators. We vary T = £/N T a a by varying A^ T . 

2) Physical problems: Increasing the temperature is 
expected to induce significant changes in the structure of 
the hadrons (see, e.g. Q for reviews). Two "extreme" 
pictures are frequently used for the intermediate and the 
high T regimes, respectively: the weakly interacting me- 
son gas, where we expect the mesons to become effective 
resonance modes with a small mass shift and width due 
to the interaction; and the quark-gluon plasma (QGP), 
where the mesons should eventually disappear and at 
very high T perturbative effects should dominate. These 
genuine temperature effects should be reflected in the low 
energy structure in the mesonic channels. But this struc- 
ture cannot be observed directly, due to the inherently 
coarse energy resolution X/l r = T. Our strategy is the 
following: we fix at T — a mesonic source which gives 
a large (almost 100%) projection onto the ground state. 
Then we use this source to determine the changes induced 
by the temperature on the ground state. For T > T c we 
do not have a good justification to use that source as 
representative of the meson but we assume that it still 
projects onto the dominant low energy structure in the 
spectral function. This is a reasonable procedure if the 
mesons interact weakly with other hadron-like modes in 
the thermal bath and the changes in the correlators are 
small. Large changes will signal the breakdown of this 
weakly interacting gas picture and there we shall try to 
compare our observations with the QGP picture. 

On a periodic lattice the contribution of a pole in 
the mesonic spectral function to the i-propagator is 
cosh(m(t — N T /2)) (this m is therefore called "pole- 
mass" ) . A broad structure or the admixture of excited 
states leads to a superposition of such terms. Fitting a 
given i-propagator by cosh(m(i)(i — N T /2)) at pairs of 
points t, t + 1 defines an "effective mass" m(t) which is 
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constant if the spectral function has only one, narrow 
peak. We shall simply speak about m(t) as "mass": it 
connects directly to the (pole) mass of the mesons below 
T c , while above T c it will help analyze the dominant low 
energy structure in the frame of our strategy above. By 
contrast, we shall speak of "screening mass" (mS *) when 
extracted from spatial propagators, to/ ct ) is different from 
the T > mass (the propagation in the space directions 
represents a T = problem with finite size effects). 

We use lattices of 12 3 x N T with N T = 72,20,16 
and 12 at = 5.68, lG = 4. We find @ T c at 
N T slightly above 18, which fixes for the above lattices 
T ~ 0, 0.93T C , 1.15T C and 1.5T C and a T ~ 0.044 fm = (4.5 
GeV) -1 . We present two sets of results: Set- A represents 
a prospective study of the T— dependence of the temporal 
propagators, calibrated with £ = 5.9(5) (apparently over- 
estimated); Set-B represents a more precise analysis of 
the T dependence of the temporal and spatial correlators 
for 3 quark masses, with also a more precise calibration 
£ ~ 5.3(1). The various parameters are given in the Ta- 
ble. Details will be given in a forthcoming paper fij. We 
use periodic (anti-periodic) boundary conditions in the 
spatial (temporal) directions and gauge-fix to Coulomb 
gauge. We investigate correlators of the form: 

G M (x,t) = Ez, yi ,y 2 W l(yi) W 2(v2) X 

(Tr [S(yx,0; z, t) lAn5 S\y 2 ,0; z + x, t)^ M ] ) (4) 
wi,2(y) ~ S(y) (point); w lt2 (y) ~ exp(-ay p ) (exp.) (5) 

with S the quark propagator and jm = {75,7ij 1,717s} 
for M — {Ps,V,S,A} (pseudoscalar, vector, scalar and 
axial- vector, respectively). We use point and smeared 
(exponential) quark sources and point sink. We fix the 
exponential (exp) source taking the parameters a, phi (|5|) 
from the observed dependence on x of the temporal Ps 
correlator with point-point source at T ~ (see Table). 
The results of a variational analysis using point-point, 
point-exp and exp-exp sources indicate that the latter 
ansatz projects practically entirely onto the ground state 
at T = 0. This is well seen from the effective mass in Fig. 
|]. Therefore we use throughout the exp-exp sources in the 
Table, according to our strategy for defining hadron op- 
erators at high T M. All masses are given in units a" 1 , 
i.e. we plot mass x a T . Errors are statistical only. 

a) Effective masses. In Fig. || we show the effec- 
tive mass m(t) of the Ps and V time-propagators at 
T ~ 0.93T C , 1.15T C and 1.5T C . The similarity between the 
Set- A and -B data indicates that calibration uncertain- 
ties are unimportant. We notice practically no change 
from T = (Fig. |) to 0.93T C , while above T c clear 
changes develop: m(t) depends strongly on t, it increases 
significantly, and the Ps and V reverse their positions. 
Because of the large changes at high T we compare here 
with the unbound quark picture of the high T regime 
of QGP. For this we calculate mesonic correlators using, 
with the same source, free quark propagators Sq instead 



of S in (g), with 7f = C = 5.9(5.3) for Set-A(Set-B) and, 
illustratively, too = m q a a — m q a T ^ = 0.1. We did not 
attempt a quantitative comparison at present but simi- 
larities are apparent - see Fig. ^. Generally above T c , 
m(t) shows stronger ^-dependence, which means that the 
spectral function selected by the exp-exp source no longer 
has just one, narrow contribution, as for T < T c . 

b) Wave functions. For the Set-B we have also an- 
alyzed the "wave functions", i.e. the behaviour of the 
temporal correlators with the quark-antiquark distance 
at the sink. In Fig. |3| we compare the Ps wave functions 
normalized at x = 0, Gp s (x,t)/Gp s (0,t), at several t 
for T ~ 0.93T C (which is very similar to T ~ 0) and 
T ~ 1.5T C at our lightest quark mass (n a = 0.086) and for 
the free quark case (m q a a = 0.1, = £ = 5.3). Our exp- 
exp source appears somewhat too broad at T ~ 0.93T C : 
the quarks go nearer each other while propagating in t. 
Interestingly enough, this is also the case at T ~ 1.5T C : 
the spatial distribution shrinks and stabilizes, indicating 
that even at this high temperature there is a tendency for 
quark and anti-quark to stay together. This is in clear 
contrast to the free quark case which never shows such 
a behaviour regardless of the source (in Fig. || we use 
the exp-exp source and too = 0.1; for heavy free quarks 
we expect a "wave function" similar to the source at all 
t) . Hence the only effect of the temperature on the wave 
function is to make it slightly broader. The same holds 
also for the other mesons at all quark masses. 

c) Masses and screening masses. On the Set-B data 
we fit the Ps, V, S and A time-correlators to single hy- 
perbolic functions at the largest 3 t- values ||. Above T c 
we assume that these masses characterize the dominant 
low energy structure selected by our source (e.g., the pu- 
tative, unstable states described by the wave functions 
discussed above). We also extract screening masses from 
spatial correlators at the largest 3 x. Since the spatial 
physical distance is large we use point-point source for 
all T (a gauge invariant extended source leads to similar 
results). The results for to and my*' at k ct = 0.086 are 
shown in Fig. || (to, and to a smaller extent also m^ 
may be overestimated). We extrapolate to and m^ in 
1/k = 2mo + 8 to the chiral limit from the 3 quark masses 
analyzed fiofl . Up to T c screening masses and masses re- 
main similar, while above T c the former become much 
larger than the latter, both at finite quark mass and in 
the chiral limit - see Fig. ||| (ll]] . In the high T regime of 
QGP one expects m {cr) ~ 2tt/N t - 0.4(0.5) in units of 
a^ 1 for T ~ 1.15T C (1.5T C ), respectively, to be compared 
with the values in Fig. |^ of ~ 0.3(0.4). We introduce: 

m ( CT ) _ m 
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as a phcnomcnological parameter to succinctly quantify 
this behaviour. Since at high T the quarks are expected 
to exhibit an effective mass m e ^ ~ gT/y6 J12j, R ~ 
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1 — 0.26c/ and thus can serve as indicator for how near we 
are to the hight T, perturbative regime. See Fig. §. 

In conclusion, in this quenched QCD analysis the 
changes of the meson properties with temperature ap- 
pear to be small below T c , while above T c they become 
important and rapid, but not abrupt. Here we observe 
apparently opposing features: On the one hand, the be- 
haviour of the propagators, in particular the change in 
the ordering of the mass splittings could be accounted for 
by free quark propagation [ p"3| in the mesonic channels 
above T Cl which would also explain the variation of m(t). 
On the other hand, the behaviour of the wave functions 
obtained from the 4-point correlators suggests that there 
can be low energy excitations in the mesonic channels 
above T c , remnants of the mesons below T c . They would 
be characterized by a mass giving the location of the cor- 
responding bump in the spectral function. The variation 
of m(t) with t and with the source would then indicate a 
resonance width increasing with T, although it may sim- 
ply reflect the uncertainty in our treatment of the low 
energy states H] . Remember that our source is not cho- 
sen arbitrarily but such as to reproduce a "pure" meson 
source at T = 0; at high T, however, it may allow admix- 
ture of other contributions, and m becomes increasingly 
ambiguous. We see chiral symmetry restoration above 
T c both in the masses and in the screening masses, with 
the latter increasing faster than the former and remain- 
ing below the free gas limit at T ~ 1.5T C . The exact 
amount of splitting among the channels and the precise 
ratio between m and may, however, be affected also 
by uncertainties in our £ calibration. Finally, note that 
this is a quenched simulation, with incomplete dynamics. 

A possible physical picture is this: Mesonic excitations 
subsist above T c (up to at least 1.5T C ) as unstable modes 
(resonances), in interaction with unbound quarks and 
gluons. Our results are consistent with this, but there 
may be also other possibilities (cf M,E6|, cf || and ref- 
erences therein). E.g., in a study of meson propagators 
including dynamical quarks but without wavefunction in- 
formation |l6f| , one found masses and (spatial) screening 
masses oc T above T c and indication for QGP with "de- 
confined, but strongly interacting quarks and gluons" . 
The complex, non-perturbative structure of QGP (al- 
ready indicated by equation of state studies up to far 
above T c |l7| ) is also confirmed by our analysis of gen- 
eral mesonic correlators. From our study however, the 
detailed low energy structure of the mesonic channels ap- 
pears to present further interesting, yet unsolved aspects. 

Further work is needed to remove the uncertainties still 
affecting our analysis. This concerns particularly the £ 
calibration and the question of the definition of hadron 
operators at high T, which appear to have been the major 
deficiencies, besides the smaller lattices, affecting earlier 
results [ pUfl . We shall also try to extract information di- 
rectly about the spectral functions fl9| . 
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FIG. 1. Effective pseudoscalar mass m(t) in units of a T 1 
vs t, for various sources at T ~ (Set-yl). 




FIG. 2. From left to right, effective mass m(t) in units of 
a^ 1 at T ~ 0.93, 1.15 and 1.5T C (open symbols) vs t. Also 
shown are the effective masses from the same correlators cal- 
culated using free quarks. A: Set-A, B: Set-B, K a = 0.086. 
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FIG. 3. Ps wave functions (Set-B, k ct = 0.086, exp-exp 
source) normalized at r = vs quark separation r at t = 2, 4 
and 6, using full and free propagators. Also plotted is the 
initial distribution of separations as given by the source, 
/ d s yw(y + r)w{y). T ~ 0.93T C (left) and T ~ 1.5T C (right). 
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FIG. 4. Temperature dependence of the mass m (open 
symbols) and screening mass (full symbols) in units a^ 1 , 
for Set-B, K a = 0.086 (upper plot) and in the chiral limit 
(lower plot). The vertical gray lines indicate T c . The data 
correspond to T ~ 0, 0.93T C , 1.15T C and 1.5T C . 
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FIG. 5. Temperature dependence of R, eq. (|6|). The 
vertical gray line indicates T c . The data correspond to 
T~0, 0.93T C , 1.15T C and 1.5T C . 
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20 


0.068 


5.4 


0.109(1) 


0.132(2) 


0.442 


1.298 
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60 


0.081 


4.05 


0.178(1) 


0.196(1) 


0.379 


1.289 
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0.084 


3.89 


0.149(1) 


0.175(1) 


0.380 


1.277 




60 


0.086 


3.78 


0.134(1) 


0.164(1) 


0.380 


1.263 



TABLE I. Simulation parameters used at every T and me- 
son masses at T ~ (in units a^ 1 ). The source parameters 
a, p eq. (^) are extracted from the T ~ wave function. 
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